Attosecond nanoplasmonic streaking of localized fields near metal nanospheres 
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Collective electron dynamics in plasmonic nanosystems can unfold on timescales in the attosec- 
ond regime and the direct measurements of plasmonic near- field oscillations is highly desirable. We 
report on numerical studies on the application of attosecond nanoplasmonic streaking spectroscopy 
to the measurement of collective electron dynamics in isolated Au nanospheres. The plasmonic field 
oscillations are induced by a few-cycle NIR driving field and are mapped by the energy of photoemit- 
ted electrons using a synchronized, time-delayed attosecond XUV pulse. By a detailed analysis of 
the amplitudes and phase shifts, we identify the different regimes of nanoplasmonic streaking and 
study the dependence on particle size, XUV photoelectron energy and emission position. The simu- 
lations indicate that the near-fields around the nanoparticles can be spatio-temporally reconstructed 
and may give detailed insight into the build-up and decay of collective electron motion. 

PACS numbers: 73.20.Mf, 78.47.J- 



Nanoplasmonics has rapidly evolved and numerous 
techniques have been developed to study the effects of 
plasmonic field enhancement^, but so far the direct, 
time-resolved measurement of plasmonic near-fields has 
not been achieved. The fastest dynamics in plasmonic 
nanosystems can take place on timescales down to 100 
attoseconds as determined from the inverse bandwidth 
of plasmonic resonance spectra. Attosecond metrol- 
ogy has provided valuable tools for measurements of ul- 
trafast electron dynamics in atomJ^^, molecules^ and 
surfaces'". One of the most successful techniques is at- 
tosecond streaking spectroscop}^^, employing photoe- 
mission of electrons by an attosecond XUV pulse syn- 
chronized to a strong optical field. Recording the elec- 
tron kinetic energy spectra (by e.g. time-of-flight (TOF) 
spectroscopy) as a function of the delay between the two 
pulses allows for the reconstruction of the laser fields and 
the electron emission dynamics. The technique is conse- 
quently a promising candidate for the real-time observa- 
tion of collective electron motion in nanosystems. 

One of the key aspects in traditional attosecond streak- 
ing on e.g. gas phase atomic targets is the spatial homo- 
geneity of the driving laser field. In contrast, noble metal 
nanoparticles exhibit strongly enhanced, but highly lo- 
calized optical fields, such that the assumption of spatial 
homogeneity is no longer valid. When applying attosec- 
ond streaking spectroscopy to nanoparticles, the spatial 
decay of the near-field into free space will govern the 
streaking process. Nanoplasmonic streaking was origi- 
nally proposed for the instantaneous electric field probing 
regim^^. A study on integrated streaking spectroscopy 
on nanostructured antennas showed, that a reconstruc- 
tion of the relatively homogeneous field in the antenna 
gap is possible if photoemission is limited to this regioiP^. 
Here the electron acceleration mostly takes place in the 
ponderomotive regime. 

We extend this previous work to spherical Au nanopar- 
ticles of different sizes and exploring the transition 
from the instantaneous to the ponderomotive streaking 
regime. By taking into account the full spatial and tem- 
poral dependence of the plasmonic near-field, we analyze 



the streaking processes as a function of the nanoparticle 
size and laser parameters in detail. Spherical particles 
allow the use of Mie theory and were chosen for compu- 
tational efficiency. The conclusions drawn in our study 
are, however, not limited to spheres, as the spatial decay 
of the near- field for a variety of systems shows similar, ex- 
ponential behavioi'^. Employing isolated nanoparticles 




FIG. 1: (color online) Scheme for attosecond nanoplasmonic 
streaking of Au nanospheres employing few-cycle NIR pump 
pulses, attosecond XUV probe pulses and the measurement 
of the kinetic energies of photoemitted electrons using a time- 
of-flight (TOF) spectrometer. The inset shows the field en- 
hancement as calculated by Mie theory for 2/ = 0. Due to 
symmetry of the plasmonic field, the point of electron emis- 
sion is uniquely defined by the angle a. 

omits the requirement for spatially resolved experiments 
such as photoelectron emission microscopy, which is ex- 
perimentally challenging in conjunction with attosecond 
light sourceJI^Ell. Such isolated nanoparticles of nearly 
equal size and shape may be prepared by chemical syn- 
thesis in large quantities'^. The nanoparticles can be 
supplied in sufficient density t o the laser focus by employ- 
ing either aerodynamic lensed^^^, optical trapping'^ or 
laser ablation'^. From the experimental point of view, 
there are two major advantages of such an approach. 
First, the contrast in the spectra is expected to be higher. 
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as this approach does not include a substrate. Further- 
more, as the particles are replaced for each laser shot, 
the damage threshold of the Au nanosphere^^ does not 
play a role and the non-linear plasmonic regime can be 
studied as well. 

The scheme for attosecond nanoplasmonic streaking of 
fields near a metal nanosphere is depicted in figure [l] 
A few-cycle NIR laser pulse excites the plasmonic oscil- 
lations and a time-delayed attosecond XUV laser pulse 
(with a photon energy on the order of 100 eV) photoemits 
electrons from the nanosphere. The photoelectrons un- 
der consideration here are emitted from the Fermi edge 
of the conduction band. Furthermore they are assumed 
to not have scattered with other electrons or the lattice, 
which is justified if emission takes place close to the sur- 
face as the electron's mean free path is on the order of a 
few monolayers'^. Accordingly the initial electron kinetic 
energy is given by Ei^i^ = hvxuv — Wf for the XUV fre- 
quency Vxuv and the work function of Au of Wf = 5.1 eV. 
While traveling to the detector, the electrons are accel- 
erated by the local electric field near the metal nanopar- 
ticle resulting in a kinetic energy shift as a function of 
delay time that is measured with a TOF spectrometer. 
In our simulations the final drift velocity of a photoelec- 
tron (emitted at time te with an initial velocity vq) was 
obtained by integrating the electron's classical equation 
of motion 
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where e is the elementary charge and m is the electron 
mass. The integration was performed with an explicit 
Runge-Kutta algorithm (4th order) using a time step of 
10 as. The plasmonic response and the corresponding 
near-field E(f^t) is calculated using Mie theory; Us- 
ing this approach, the enhancement and phase shift of 
the electric field outside the sphere can be calculated by 
a normal mode expansion of the sphere's electromagnetic 
response. The incident, linearly polarized laser field is 



E{t) = Eoexp 



cos {ut + 0) e^;. 
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with r = 2.1 fs, corresponding to a full- width at half- 
maximum of 5 fs, an angular frequency uj corresponding 
to a center wavelength of A = 720 nm, and a carrier- 
envelope phase ^ = TT. The Coulomb field of the charged 
nanoparticle after electron emission is orders of magni- 
tude smaller than the plasmonic near-field and thus ne- 
glected. 

Using the appropriate wavelength dependent dielectric 
function for Au fitted to experimental data, Mie theory 
predicts the resonance of Au spheres to be well below 
600 nm for diameters up to 100 nm (see Suppl. Inf.). In 
turn, for these diameters the system is driven off reso- 
nance with the implication that the field enhancement as 
well as the phase of the enhanced near-field do not vary 
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FIG. 2: (color online) a) Phase shift /\(j)s and b) maximum 
velocity shift as a function of electron energy Ekin and Au 
sphere diameter d obtained for a laser intensity of 10^^ W/cm^ 
and a = 0°. 



significantly over the entire NIR laser spectrum. Under 
this condition, the calculations can be simplified and the 
entire laser spectrum treated with a single set of Mie coef- 
ficients for the center wavelength. The phase shift A(/)p of 
the near-field with respect to the driving field was taken 
into account in the calculations, but plays only a minor 
role for small Au spheres. The field enhancement factor 
(J = Ep/Eq for Au spheres is predicted to be between 3.4 
and 4 for spheres with diameters d between 10 nm and 
100 nm, respectively. Throughout this study we use pa- 
rameters for which the maximum velocity shift of a pho- 
toelectron is small compared to its initial velocity, so that 
recollision and rescattering do not have to be taken into 
account. Electrons resulting from multi-photon emission 
and acceleration in the NIR field are neglected since for 
the chosen NIR intensity they are energetically well sep- 
arated from the XUV emitted electrons^^. 

First, we limit the analysis to electrons liberated from 
a = 0° (see fig. [T]) and traveling in x-direction corre- 
sponding to the emission along the polarization axis to- 
wards the TOF spectrometer. Significant insight can be 
gained from the phase shift between the streaking curve 
and the plasmonic near-field A^^. The phase shift serves 
as a measure how the streaking spectrum is generated 
and is in principle dependent on all experimental pa- 
rameters (namely sphere diameter, electron energy, en- 
hancement factor, duration of the near- field). Figure [2] 
a) shows A(j)s as a function of the sphere diameter d and 
electron energy E^in for r = 2.1 fs and an intensity of 
lO^^W/cm^ of the NIR field. The size dependence of 
A(j)s originates from the spatial extension of the near- 
fields into the vacuum, which is decaying with approxi- 
mately 1/r^, with r being the radial distance from the 
origin. For moderate electron energies and large spheres, 
the phase shift is approaching 7r/2. 

In this regime the electron velocity is small with re- 
spect to the field decay length, resulting in a pondero- 
motive acceleration in the plasmonic field. Decreasing 
the particle size at a fixed initial kinetic energy leads 
to a decrease of A^^. Now the regime of instantaneous 
near-field probing is approached, as the electrons leave 
the near-field faster and the acceleration is thus also con- 
fined in time. Going to even smaller particles, A<ps even- 
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tually increases again, because here the net acceleration 
by the near-field is smaller than the ponderomotive accel- 
eration in the laser field. The minimum achievable phase 
shift is mainly determined by the enhancement factor <j, 
which defines the relative strength of the two contribu- 
tions. The maximum velocity shift Avmax of the streak- 
ing curve is shown in Figure |2] b). For the moderate 
near-field intensities used here, the velocity shift increases 
linearly with the electric near-field amplitude. Conse- 
quently, Avmax generally follows the field enhancement 
factor peaking for a sphere diameter of about 170 nm. 
For small sizes and/or large electron energies, Avmax ap- 
proaches the value expected from streaking in the ho- 
mogenous laser field of 2 x 10~^m/s. 

For the spherical Au particles under consideration with 
relatively low enhancement factors, we can conclude that 
a direct probing of the electric field, as was discussed 
in Stockman et al}^ is not achievable. On the other 
hand, the enhancement factor is large enough to give 
sufficient contrast in streaking spectrograms especially 
in the ponderomotive regime. 
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FIG. 3: (color online) Streaking traces for emission positions 
a = — 90° for 10 nm (a) and 100 nm spheres (b) and the 
respective electric fields acting on the photoelectrons (emitted 
at te = fs) as a function of time after photoemission, (c) and 
(d). The blue dashed lines show the plasmonic near field. 

Another important aspect arising from the spatial in- 
homogeneity of the plasmonic fields is the dependence 
of nanoplasmonic streaking on initial conditions such as 
emission position and launch angle. This is especially im- 
portant, because XUV radiation can liberate electrons at 
any position on the sphere. For the following discussion 
the XUV photon energy is kept fixed at hvxuv = 105 eV, 
which can routinely be produced via high-harmonic gen- 
eration in Ne^^. Figures [3^) and b) show the infiuence 
of the emission position a on the streaking traces for 10 
nm and 100 nm Au spheres, respectively. In general, for 
larger values of a (representing emission from a ring) the 
phase shift A^^ increases while the streaking amplitude 
decreases. This can be explained by the field distribu- 
tion around the nanosphere (see figure [T]). Note that the 
particle sizes considered here are much smaller than the 



wavelength, such that the dipolar mode is dominating 
and all spheres show qualitatively similar field distribu- 
tions. 

Neglecting additional acceleration, the 100 eV photo- 
electrons travel about 6 nm per femtosecond. Electrons 
emitted on-axis {a = 0°) will show the largest energy 
shift and smallest phase shift A(/)s, as they are acceler- 
ated by the strongest near-field. Photoelectrons emit- 
ted at the sides of the sphere {a ^ 90°) initially experi- 
ence weak or even opposite electric fields (with respect 
to a = 0°) after liberation. The electric fields acting on 
electrons that have left the sphere from different emis- 
sion positions and travel to the detector will eventually 
converge. Figure (sj^) and d) show the effective, time- 
dependent electric fields acting on electrons for 10 nm 
and 100 nm Au spheres, respectively. For the 10 nm 
sphere the fields of all emission positions have become 
equal after less than one half-cycle of the plasmonic field 
oscillation. For larger spheres this situation changes sig- 
nificantly. The electron travel time through the near-field 
approximately scales with the sphere radius. Electrons 
emitted from the side of e.g. a 100 nm sphere do not 
reach the regions of strong field enhancement before the 
near-field oscillations vanish (owing to the short pulse 
duration discussed here). Accordingly, these electrons 
do not show noticeable net acceleration. 

The complex dependence of the electron energy shift on 
the initial emission conditions will determine nanoplas- 
monic streaking spectrograms recorded under the de- 
scribed geometry. A distortion of the streaking spec- 
trograms especially for larger sizes is expected from the 
analysis above. To study these effects in more detail, 
a large number of electron trajectories with randomized 
initial conditions were integrated for each sphere size. 
The initial conditions accounted for the probabilities of 
the emission position, launch angle (within a 44° cone, 
corresponding to the acceptance angle of the TOF spec- 
trometer) , kinetic energy within the XUV spectral band- 
width (5 eV) and XUV pulse duration (200 as). Trans- 
mission of XUV through the nanosphere and emission 
from the back side was included in the calculations (see 
Suppl. Inf.). The trajectory calculations were performed 
for the few-cycle NIR pulses described above (r = 2.1 fs, 
/ = 10i2^/cm2,^ = 7r). 

Figures |4^) and b) show streaking spectrograms com- 
posed of 1.5x10^ electron trajectories for 10 nm and 
100 nm spheres, respectively. The two streaking spec- 
trograms are vastly different. The amplitude and phase 
shifts for different emission positions described above are 
responsible for this behavior. The typical streaking pat- 
tern is almost completely washed out in |4]3) (100 nm), 
while it can still be observed in |4^) (10 nm). Electrons 
emitted from a 100 nm sphere at angles a larger than 
60° are hardly accelerated. Since the statistical weight 
of each contribution to the spectrogram is approximately 
proportional to the respective area on the sphere, around 
50% of the photoemitted electrons exhibit minimal en- 
ergy shifts. 
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FIG. 4: (color online) Simulated streaking spectrograms for 
10 nm (a) and 100 nm (b) Au spheres. The red line shows the 
contribution of the a = 0° trajectories. 

A streaking spectrogram recorded for a rare gas (e.g. 
Ne)^ under the same experimental conditions can serve 
as a reference permitting the precise reconstruction of the 
driving NIR field and XUV pulse parameter^. Such a 
reconstruction is also desirable for the temporal and spa- 
tial properties of the plasmonic field. As the equation 
of motion of the electrons in its general form (eq. [T]) is 
not analytically integrable, numerical simulations as pre- 
sented in Fig. [4] have to be performed for the retrieval of 
the spatio-temporal near-fields of the nanospheres. The 
reconstruction could employ a feedback algorithm, which 
modifies the simulation parameters until best agreement 
with the measured spectrogram has been achieved. 

Under the assumption of a constant average electron 
velocity and approximating the spatial near-field decay 
by an exponential function, an analytical integration of 
the equation of motion for a = 0° trajectories is possible 
(see Suppl. Inf.). This allows for a computational less 
expensive retrieval. A major difficulty, however, lies in 
extracting this specific contribution from the streaking 
spectrogram. As can be seen in figure |4j the streaking 
waveform for a = 0° (red line) can be identified more 



easily for the 100 nm sphere and enables an analytic re- 
construction. 

In order to test the experimental implementation of 
attosecond nanoplasmonic streaking, we estimated the 
expected electron count rates using table-top attosecond 
XUV light sources and typical nanoparticle target den- 
sities provided by aerosol evaporation techniques^. Al- 
though we have only considered single particles in our 
simulations, an ensemble of nanoparticles may be utilized 
in an experiment for a target that provides sufficient spa- 
tial separation between individual particles. For an XUV 
pulse with 10^ photons, a target density of 10^ cm~^ and 
100 nm Au particles we obtain a count rate of 1 s~^ using 
a 10 kHz laser system (see Suppl. Inf.). The count rate 
can be increased by one-three orders of magnitude by 
employing optical trapping methods^^ or laser ablation 
targets^. Alternatively, an XUV light source with very 
high flux, such as the one expected to emerge from the 
European Light Infrastructure (ELI)^^, could be utilized. 

In conclusion, we analyzed and discussed the attosec- 
ond streaking of plasmonic near-fields around isolated Au 
nanospheres of different diameters. In contrast to con- 
ventional attosecond streaking, the electron acceleration 
and deceleration has a strong spatial dependence and sig- 
nificant differences are identified in our simulations for 10 
nm and 100 nm particles. These differences are caused by 
retardation effects between electron trajectories emitted 
from different positions on the nanosphere. The spatio- 
temporal information that can be retrieved from recorded 
streaking spectrograms could greatly deepen our under- 
standing of how plasmons form and how they decay. In 
contrast to studies on nanostructured surfaces, where the 
field strength must remain below the damage threshold 
of the material, a beam of isolated nanoparticles would 
allow studies in the non-linear regime. Furthermore, the 
approach can be extended to other material groups and 
be used to e.g. study - in real-time - ultrafast field- 
induced phase transitions, su ch as the metallization of 
semiconductors and dielectric d^^ l ^^l 
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